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We give the nonperturbative phase diagram of the four-dimensional hot electroweak phase transition. A sys- 
tematic extrapolation a — > is done. Our results show that the finite temperature SU(2)-Higgs phase transition is 
of first order for Higgs-boson masses ran < 66.5±1.4 GeV. The full four-dimensional result agrees completely with 
that of the dimensional reduction approximation. This fact is of particular importance, because it indicates that 
the fermionic sector of the Standard Model (SM) can be included perturbatively. We obtain that the Higgs-boson 
endpoint mass in the SM is 72.4 ± 1.7 GeV. Taking into account the LEP Higgs-boson mass lower bound excludes 
any electroweak phase transition in the SM. 



The observed baryon asymmetry is finally de- 
termined at the electroweak phase transition 
(EWPT) jjj. The perturbative approach breaks 
down for the physically allowed Higgs-boson 
masses (e.g. TOh > 70 GeV) M. Since merely the 
bosonic sector is responsible for the bad perturba- 
tive features (due to infrared problems) the simu- 
lations are done without the inclusion of fermions 
on four-dimensional lattices ||, [Q. Another ap- 
proach is the simulations of the reduced model 
in three-dimensions The comparison of 

the results obtained by the two techniques is not 
only a useful cross-check on the perturbative re- 
duction procedure but also a necessity because 
the fermions must be included perturbatively. 

The line separating the symmetric and broken 
phases on the m# — T c plane has an endpoint, 
rriH,c- In four dimension at mn ~ 80 GeV the 
EWPT turned out to be extremely weak, even 
consistent with the no phase transition scenario 
on the 1.5-er level 0). Three-dimensional results 
show that for tub > 95 GeV no first order phase 
transition exists ^] and more specifically that the 
endpoint is mH, c « 67 GeV §. In this paper 
we present the analysis of the endpoint on four 
dimensional lattices in the continuum limit and 
transform the result to the full SM. 

We will use different spacings in temporal (at) 
and spatial (a s ) directions. The asymmetry of the 
lattice spacings is given by the asymmetry factor 



£ = a s /a t . The action reads 

s\um = 
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where U x ^ denotes the SU(2) gauge link variable, 
U sp and Utp the path-ordered product of the four 
U x ,n around a space-space or space-time plaque- 
tte, respectively; ip x stands for the Higgs field. 
The anisotropies 7I = Pt/Ps and 7^ = Kt/n s 
are functions of the asymmetry £. These func- 
tions have been determined perturbatively |l0| 
and non-perturbatively [|ll| . In this paper we use 
the asymmetry parameter £ = 4.052, which gives 
j K = 4 and 73 = 3.919. 

We have performed our simulations on finer 
and finer lattices, moving along the lines of con- 
stant physics (LCP). In our case there are three 
bare parameters (k,/3, A). The bare parameters 
are chosen in a way that the zero temperature 
renormalized gauge coupling gn is held constant 
and the mass ratio for the Higgs- and W-bosons 
Rhw = in-H I mw corresponds to the Higgs mass 
at the endpoint of first order phase transitions: 
Rhw,c- These two conditions determine a LCP as 
a one-dimensional subspace in the original space 
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of bare parameters. The position on the LCP 
gives the lattice spacing a. As the lattice spacing 
decreases Rhw,c — > Rnw^ont.- 

Since our theory is a bosonic one we assumed 
that the corrections are quadratic in the lattice 
spacings; thus an a 2 fit has been performed for 
Rhw,c to determine its continuum value. 

We have carried out T ^ simulations on 
L t = 2,3,4,5 lattices and tuned k to the tran- 
sition point. This condition fixes the lattice spac- 
ings: a t — a s /^ = l/(T c L t ) in terms of the tran- 
sition temperature T c in physical units. The third 
parameter A, finally specifying the physical Higgs 
mass in lattice units, has been chosen in a way 
that the transition corresponds to the endpoint 
of the first order phase transition subspace. 

In this paper V = Lt - l? s type four-dimensional 
lattices are used. For each L t we had 8 different 
lattices, each of them had approximately twice 
as large lattice-volume as the previous one. The 
smallest lattice was V = 2 • 5 3 and the largest one 
was V = 5 ■ 50 3 . Reweighting was used to obtain 
information in the vicinity of a simulation point. 

The determination of the endpoint of the finite 
temperature EWPT is done by the use of the Lee- 
Yang zeros of the partition function Z . Near 
the first order phase transition point the parti- 
tion function reads Z = Z s + Z\> cx exp(— Vf s ) + 
exp(-Vfb) , where the indices s(b) refer to the 
symmetric (broken) phase and / stands for the 
free-energy densities. Near the phase transition 
point we also have /(, = f s + a(n — n c ) , which 
shows that for complex k values Z vanishes at 
Im(n) = 7r • (n — 1/2) /{Va) for integer n. In case 
a first order phase transition is present, these Lee- 
Yang zeros move to the real axis as the volume 
goes to infinity. In case a phase transition is ab- 
sent the Lee- Yang zeros stay away from the real 
k axis. Denoting kq the lowest zero of Z, i.e. 
the position of the zero closest zero to the real 
axis, one expects in the vicinity of the endpoint 
the scaling law Im(« ) = ci(L t , \)V + C2{L t ,X). 
In order to pin down the endpoint we are look- 
ing for a A value for which ci vanish. In practice 
we analytically continue Z to complex values of 
k by reweighting the available data. Also small 
changes in A have been done by reweighting. As 
an example, the dependence of ci on A for L t — 3 



is shown in fig. 2. To determine the critical value 
of A i.e. the largest value, where C2 = 0, we have 
performed fits linear in A to the nonnegative Ci 
values. 
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Figure 1. Dependence of C2 on A for L t = 3. 



Having determined the endpoint X C rit.(-^t) for 
each Lt we calculate the T = quantities 
(Rhw,9 2 r ) on V = (32Z t ) ■ (8i t ) ■ (6L t ) 3 lattices, 
where 32Lt belongs to the temporal extension, 
and extrapolate to the continuum limit. Having 
established the correspondence between \ C rit.{L t ) 
and Rhw > t ne L t dependence of the critical Rhw 
is easily obtained. Fig. 3 shows the dependence 
of the endpoint Rhw values on l/L 2 . A linear 
extrapolation in l/L 2 yields the continuum limit 
value of the endpoint Rhw- We obtain 66.5 ± 1.4 
GeV, which is our final result. 

Comparing our result to those of the 3d anal- 
yses m one observes complete agreement. Since 
the error bars on the endpoint determinations are 
on the few percent level, the uncertainty of the 
dimensional reduction procedure is also in this 
range. This indicates that the analogous per- 
turbative inclusion of the fermionic sector results 
also in few percent error on My. 

Based on our published data ji 11 and the 
present results we are able to draw the phase dia- 
gram of the SU(2)-Higgs model in the (T c /niH — 
Rhw) plane. This is shown in fig. 4. The con- 
tinuous line - representing the phase-boundary - 
is a quadratic fit to the data points. 

Finally, we determine what is the endpoint 
value in the full SM. Our nonperturbative anal- 
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Figure 2. Dependence of Rnw.cr, i-e. Rhw cor- 
responding to the endpoint of first order phase 
transitions on 1/Lf and extrapolation to the in- 
finite volume limit. 
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Figure 3. Phase diagram of the SU(2)-Higgs 
model in the (T c /rriH — Rhw) plane. 



ysis shows that the perturbative integration of 
the heavy modes is correct within our error 
bars. Therefore we use perturbation theory |L3| 
to transform the SU(2)-Higgs model endpoint 
value to the full SM. We obtain 72.4 ± 1.7 GeV, 
where the error includes the measured error of 
Rnw,cont., 9r and the estimated uncertainty jL4} 
due to the different definitions of the gauge cou- 
plings between this paper and |fl3f . The dominant 
error comes from the uncertainty on the position 
of the endpoint. 

The present experimental lower limit of the 
SM Higgs-boson mass is 89.8 GeV ||. Taking 
into account all errors (in particular those coming 
from integrating out the heavy fermionic modes), 



our endpoint value excludes the possibility of any 
EWPT in the SM. This also means that the SM 
baryogenesis in the early Universe is ruled out. 

For details of this analysis see . 

Simulations have been carried out on the Cray- 
T90 at HLRZ-Jiilich, on the APE-Quadrics at 
DESY-Zeuthen and on the PMS-8G PC-farm in 
Budapest. This work was supported by OTKA- 
T016240/T022929 and FKP-0128/1997. 
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